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Coupling between axial and torsional degrees of freedom often modifies the con- 
formation and expression of natural and synthetic filamentous aggregates. Recent 
studies on chiral single-walled carbon nanotubes and B-DNA reveal a reversal in 
the sign of the twist-stretch coupling at large strains. The similarity in the re- 
sponse in these two distinct supramolecular assemblies and at high strains suggests 
a fundamental, chirality dependent non-linear elastic behaviour. Here, we seek the 
link between the microscopic origin of the non-linearities and the effective twist- 
stretch coupling using energy-based theoretical frameworks and model simulations. 
Our analysis reveals a sensitive interplay between the deformation energetics and 
the sign of the coupling, highlighting robust design principles that determine both 
the sign and extent of these couplings. These design principles have been already 
exploited by Nature to dynamically engineer such couplings, and have broad im- 
plications in mechanically coupled actuation, propulsion and transport in biology 
and technology. 

Keywords: chirality, supramolecular assemblies, mechanical coupling, Elastica, 
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1. Introduction 

Structural chirality (i.e. handedness) often results in large mechanical couplings in 
filamentous supramolecular assemblies. As a classical example, twist-stretch elas- 
ticity in DNA is vital during chromatin organization, transcription regulation and 
binding specificity during DNA-complex formation (Hagerman 1988; Tsai et al. 
1978; Hogan & Austin 1987; Romberg & Lorch 1995; Mirkin 2000). The morphol- 
ogy of the helical, rotating flagella that provide thrust to several classes of bacteria 
during their chemotaxis is in part dependent on the sign and extent of such mechan- 
ical couplings (Ritao et al. 2006; Coombs et al. 2002). They are also expected to 
influence the mechanics of twisted biopolymer aggregates such as coiled-coils (Offer 
et al. 2002), amyloid fibrils (Jimenez et al. 2002), and network forming fibers such 
as actin, fibrin (Weisel et al. 1987; Medved et al. 1990) and collagen (Bozec et al. 
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2007). Similar chiral supermolccular filaments have also been synthesized, in partic- 
ular at the nanoscale. In applications which rely on physico-mechanical couplings, 
engineering the coupled response in structurally robust helical nanowires (Kondo & 
Takayanagi 2000), nanobelt coils (Pan et al. 2001), chiral nanotubes, and their as- 
semblies offers an elegant route for fabrication of nanoscale motors, oscillators and 
switches by eliminating the need for externally actuated deformation (Craighead 
2000; Baughman et al. 1999; Baughman et al. 2002; Fenimore et al. 2003; Papadakis 
et al. 2004; Liang & Upmanyu 2005; Liang & Upmanyu 2006; Cohen-Karni et al. 
2006). 

Twist-stretch coupling in chiral, elastically stiff filaments is usually described 
by linearized Elastica (Love 1944). Well below the relevant persistence lengths, the 
effect of thermal fluctuations is negligible and the free energy per unit length of a 
straight filament subject to an axial force F can be expressed in terms of the strain 
components, 



U 



SU b (R) + l -k,e 2 + ^k^ 2 + k n <t>e+... 



Fe. (1.1) 



The first term is the contribution, if any, due to the change in radius of the filament. 
For a hollow tube of radius R, this contribution can be expressed as a change in 
the bending energy of the rollcd-up sheet due to its mean curvature k — l/R, i.e. 
Ub = \kbjR 2 with kb the bending rigidity of the rolled up sheet. The constants k e 
and k<p are the extensional and torsional rigidities of the rod, while kn couples the 
twist per unit length (f> to axial strain e. The magnitude of kn is a measure of the 
interplay between structural chirality and the microscopic deformation mode. 

Minimizing the free energy at a fixed applied strain (e or <j>) yields the equi- 
librium conjugate strain, and therefore the effective coupling. As an example, for 
negligible change in the filament radius, the coupling reduces to 

fell , k U (j)-F 

and e eq = , (1.2) 



where (f> eq and e eq are the equilibrium strains induced by prescribed conjugate 
strains e and <p, respectively. Similar energy-based models have been used to de- 
scribe the coupled mechanical response of B-DNA (Marko 1997; Kamian et al. 1997; 
Bustamante et al. 2000; Lionnct et al. 2006), twist-storing semiflcxible biopoly- 
mers (Moroz & Nelson 1998), chromatin fibers (Ben-Haim et al. 2001), and bacteria 
flagella (Coombs et al. 2002). 

This article is motivated by the recent studies on axially strained chiral single- 
walled carbon nanotubes (SWCNT) and B-DNA. Surprisingly, the effective cou- 
pling d(j)/de exhibits a sign reversal in both filaments, past a critical strain e c - 
under compression in chiral SWCNTs (Liang & Upmanyu 2006) and in stretched 
B-DNA (Gore et al. 2006). Evidently, the strain dependent coupled response cannot 
be described by linearized Elastica (equations 1.1 & 1.2). In an attempt to under- 
stand and quantify this behaviour, we isolate the microscopic response in each case, 
which then serves as the physical basis for developing energetic frameworks that 
capture the non-linear response. 
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2. SWCNT Response 

The atomic-scale simulations on the SWCNTs were performed using the second gen- 
eration reactive bond order (REBO) potential, which reproduces several mechanical 
properties of fullerenes, graphenes and diamond to a high degree of accuracy (Bren- 
ner et al. 2002). The results discussed here are limited to the (8,4) chiral nanotube, 
i.e. the SWCNT axis is set by the chiral vector Ch = 8ai + 4&2, where ai and &2 
are the graphite basis vectors. A Langevin bath was used to fix the temperature 
(0.1° K) and the strain was applied in increments of 0.25% by displacing the end 
atoms (shown masked in figure 1A). A conjugate gradient relaxation was performed 
between each successive strain step. During relaxation, the end atoms were allowed 
to displace along the tangential and radial directions. The reported twist (per unit 
length) is the relative rotation between the SWCNT ends averaged over the net 
rotation of the end atoms. For details, see Liang & Upmanyu (2006). 




Axial Strain s (%) 

Figure 1. Twist stretch elasticity in a chiral SWCNT: (top) Axial distribution of 
relative twist (color scale) in an axially strained chiral (8,4) SWCNT (Liang & Upmanyu 
2006). The helical arrangement of the graphene structural units is shown via bold lines, 
(bottom) The strain dependence of the coupling (red line) and one of the three unique 
bond angles (3 (gray line), both of which exhibit a reversal at e m 2%. (bottom left inset) 
Schematic illustrating the role of inter-fiber shear instability in reversing the coupling 
between axial and shear strains. The fibers connecting the two filaments can deform via 
fiber extension and inter-fiber shear (slip). In (a), the compression of the filaments proceeds 
via a combination of both, resulting in reorientation of the fibers and shear to the left. 
The shear is reversed in the case of enhanced inter-fiber with negligible fiber rerientation, 
as illustrated in (b). 



The sign of the SWCNT coupling constant k%i is positive at small strains and 
consistent with the SWCNT handedness (see figure 1), the chiral angle between the 
SWCNT axis and the basis vectors of the underlying graphitic lattice. The "twist- 
hardening" under compression forces the nanotube to untwist and then twist in 
the opposite direction. The change in its radius is negligible during this reversal, 
consistent with its relatively small Poisson's ratio (Yakobson et al. 1996; Jin & 
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Yuan 2003), v <~ 0.2. The extensional and torsional rigidities are weak functions 
of the SWCNT chirality (Vaccarini et al. 2000, Jin & Yuan 2003), ruling out the 
direct effect of changes in chirality due to the induced twist. 

A potential source for the non-linear response is the chirality dependence of 
the coupling constant, k\\ — fc^ sin 6(6> ± </>), where 8o is the chiral angle of the 
relaxed SWCNT and is a constant (Liang & Upmanyu 2006). This specific 
form satisfies the symmetry constraints, in particular the six fold symmetry of the 
underlying graphitic lattice, and leads to terms in higher powers of <p and a strain 
dependent coupling, 

fcn = sin 6#o ± 60 cos 6#o — 18(/. 2 sin66»o + .... (2.1) 

However, the effect of these higher order terms is small for the high chirality SWC- 
NTs for which the coupling has been observed (9 a ~ tt/6) as cos60 o ~ and 
the ratio e/cfi ~ 10 over the sign reversal range, effectively ruling out chirality de- 
pendence of the coupling constant as the basis for the observed strain dependent 
coupled response. 

Microscopic analyses of the strain dependence of the bond lengths and bond 
angles within the simulations reveal that the reversal coincides with that in the 
trend associated with one of the three bond angles between the carbon atoms 
on the SWCNT surface (angle /3, see figure 1). This change in the mechanism 
of strain accommodation effectively softens the nanotube surface with respect to 
shear, which modifies the otherwise affine response of the graphene network. The 
helically coiled graphene elements slide past each other and the enhanced surface 
shear counteracts the tendency to increase the SWCNT handedness. An analogous, 
two-dimensional response is illustrated schematically in figure 1 (inset, bottom left) 
for a sheet with bare elastic anisotropy which now couples the compression to its 
shear. The deformation of the anisotropic layer responsible for the coupled response, 
in particular the ratio of layer compression to intra-layer slip, controls the direction 
of the shear. 

The basis for this effect is an additional elastic anisotropy term in equation 1.1 
that we have ignored so far because it is weakly chirality dependent (Gartstcin et 
al. 2003), 

SU(9,<p,e) cx (e 2 -</> 2 )cos6((9 o + 0). (2.2) 

Cleary this term becomes important at large axial strains, with the leading contri- 
bution of the form ki 2 (f>e 2 , with k\ 2 a constant. Incorporating this term and max- 
imizing the resultant energy functional, we recover a strain dependent equilibrium 
twist and the critical axial strain associated with the twist reversal, 

efcn + e 2 fci2 fcn 

<j>eq = 7 . e ^~^Tr- ( 2 - 3 ) 

fc0 2fci2 

For the (8,4) SWCNT, fc = 280nNnm 2 , k u = lOOnNnm and e c = -1.3% (Vac- 
carini et al. 2000; Liang & Upmanyu 2006), which yields the higher order coupling 
constant, k\ 2 — — 3.9/xNnm. 

The strain dependent coupling points to an elegant design principle whereby 
the imposed strain can be used to engineer the coupling by effecting transitions 
in the microscopic deformation. Can we then fashion the structure so as to evoke 
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Figure 2. Twist stretch elasticity in a model chiral nanotube: Finite element simu- 
lations of twist-stretch elasticity in a model chiral filament formed by rolling a hexagonal 
lattice of I-shaped beam elements with varying beam geometries, as parameterized by 
£ = //(AL 2 ). 



the desired coupled response? As an answer, we study the coupled response in 
an identically structured model (8,4) chiral nanotube with simplified interactions. 
The carbon-carbon bonds are replaced by I-shaped beam elements of equal length 
(L = 0.142nm). The area of cross section is fixed at A = 0.005 nm 2 . The Young's 
modulus and Poisson's ratio of the beam element material are fixed at E = 1 GPa 
and v — 0.2. These elements can either stretch or bend in order to absorb the 
imposed axial strain; the twisting of the beams is negligible and ignored. The energy 
stored in bending and stretching scales quadratically with the deformation variables, 
e and the beam curvature k. The shape of the beams is important as it allows us 
to independently control the energetic competition between beam bending and 
stretching, the dimensionless ratio £ = I/(AL 2 ), where / is the moment of inertia 
of each beam. 

The software package ANSYS is used to perform the finite element simulations 
of the strained filament. A non-linear solver is employed to ensure that the errors 
at high axial strains are minimal. The axial strain and measurement of the induced 
twist is performed as in the SWCNT case. In addition, strain dependence of the 
deformations of the individual beams are also monitored. Simulations for varying 
values of £ indicate that the small strain coupling is consistent with the chirality, i.e. 
the twist changes sign with the axial strain (figure 2). Remarkably, the large strain 
response is sensitive to the value of £ and is quite rich. It varies from no reversal 
(£ = 0.04), to reversal under compression (£ = 0.03) and even tension (£ = 0.07). 

In order to gain a more rigorous understanding of the £ dependence of the 
coupled response, we have analysed the microscopic energetics of the individual 
beam elements as a function of the applied strain. Figures 3A (£ = 0.03) and 
3B (£ = 0.07) show plots of energies associated with bending and stretching the 
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Figure 3. Beam energetics in the finite element simulations: Strain dependence of 
the energy stored in beam bending and stretching for £ = 0.03 (A) and £ = 0.07 (B) model 
nanotubes. The individual beams a, b and c for which the strain are calculated (right) are 
shown in the inset. The double arrows indicate the direction of axial strain. For clarity, 
the strain interval over which the reversal in the coupled response occurs is shaded. 



beam elements for the entire filament as well as for the three unique beam elements 
(a, b and c) within a unit hexagonal cell on the filament surface. For the low £ 
filament, change in the coupling sign (shaded region) follows a gradual change 
in the energetic competition between bending and stretching under compression, 
with main contributors to this change being beams b and c. A similar transition is 
observed under tension for the high £ filament. The balance in the energy stored in 
beam c distinctly shifts towards stretching at about the same strain level at which 
we sec the reversal in the coupling for the entire filament (shaded region). Finally, 
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no such reversal in trends is observed for the intermediate £ (= 0.04) filament (not 
shown), where the coupling more or less remains constant at all strains (figure 2). 

The fact that this simplified model captures the essence of the SWCNT response 
indicates that the behaviour is both scale and model independent. The variations 
in the coupling are a direct reflection of the transitions in the microscopic energetic 
balance, quite like the observed surface shear instability in SWCNTs. Each tran- 
sition disrupts the trends in the microscopic distortions, changing the sign of the 
coupling. While the small strain response is controlled by the chirality, the strain 
dependence of the coupling arises due to atomic-scale transitions in the deforma- 
tion modes. The effect is captured in an Elastica-based framework by a higher 
order coupling constant fci2(£) which depends on the microscopic energetic balance 
and therefore modifies the energetic coupling between stretching and twisting the 
filament. 

3. B-DNA Response 




Axial strain £ (%) 

Figure 4. Twist stretch elasticity in B-DNA (left): The strain dependence of the 
twist-stretch coupling predicted by Eq. 2 (solid black line). The behaviour observed in 
bead assay experiments on a 8.4 kb long DNA molecule is also plotted for comparison 
(solid circles, Gore et al. 2006). (right) Strain dependence coupling constant fcn, predicted 
by our analysis (solid line) and that by the toy model employed by Gore et al. (2006). 

In B-DNA, the sign reversal of the effective twist-stretch coupling observed in ro- 
tor bead experiments by Gore et al. (2006) on sufficiently long individual molecules 
occurs past the entropic elastic regime. The bend fluctuations are frozen and the 
molecule behaves like a chiral elastic rod (Kamian et al. 1997; Bustamante et al. 
2000) . The critical strain associated with the reversal is well below that correspond- 
ing to the force at which the molecule undergoes a structural transition (Smith et 
al. 2006; Bustamante et al. 2000; Lionnet et al. 2006; Wereszczynski & Andricioaei 
2006). The microscopic basis for sign reversal appears to be different from that in 
the SWCNTs as the response may involve a substantial decrease in the DNA radius 
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owing to its large Poisson's ratio (y <~ 0.5), i.e. R w R a (l — i>e) with i? the relaxed 
effective radius (Gore et al. 2006; Lionnct et al. 2006). 

Then, when the extension is prescribed, the coupled response is modified due to 
radius dependence of torsional rigidity kj,, in addition to that of the coupling param- 
eter fen. Ignoring the direct effect of chirality change, the torsional rigidity scales 
as R 4 . For short-range elastic interactions, scaling and symmetry considerations re- 
quire that for harmonic distortions the coupling kn constant is an odd function of 
R, i.e. fcn(-R) = J2 n >o ^2n+i/R 2n+1 , where K 2n +i are constants. In general, these 
constants will also be function of the helicity of the DNA which may change with 
the change in radius. We ignore this relatively smaller effect. We further assume 
that the double-stranded DNA behaves as a chiral rod composed of a material with 
bare surface anisotropy, i.e. the surface shear associated with the twist r(cx 4> eq R) 
is non-zero for R — > oo. As comparison, note that in the case of SWCNTs, the 
leading order term must be 1/i? 3 as the coupling is purely curvature derived, yet 
the graphite sheet is isotropic (Liang & Upmanyu 2006) . This is a generalization of 
the toy model assumed by Gore et al. (2006), where the chirality-dependent small 
strain response was based on stiffer outer helical wire encompassing a compressible, 
isotropic inner core. 

Minimizing the energy dU/d<fi at constant axial strain e, we get 

^ = -^ ii( * )= -^(f) (§ + §+•••)' 

where Rq is the effective radius for the unstrained double strand (e = 0). Plugging 
in the strain dependence of the radius R — i?o(l — ve), a series expansion limited 
to 0[e] 3 yields the induced twist, 

<Me) = - ^ (k* n + k* 12 e + k* 13 e 2 + 0[ef) , (3.1) 

with fc*! = fcn(0), and k\ 2 & ^13 functions of Ro, v, K\ and K%. Using the known 
mechanical properties of B-DNA (R = 0.924 A, fc^(O) = 460pNnm, fcu(O) = 
— 90pNnm) and the value of the critical strain associated with sign reversal, e c = 
2.5%, we get K x = 1761 pNnm 2 and K 3 = -1433 pNnm 4 (Bustamantc et al. 2000; 
Gore et al. 2006). The predicted equilibrium twist as a function of the stretch is in 
excellent agreement with the experimental data, plotted in figure 4. However, the 
agreement rapidly breaks down past the critical strain, possibly due to increasing 
anharmonicity in the deformation (Wiggins et al. 2006) and/or an increase in the 
chirality dependence of the coupling constant. 

The analysis reveals that the magnitude of coupling constant (also plotted) de- 
creases more rapid than that predicted by Gore et al. (2006). In effect, stretched 
B-DNA always decreases its chiral character. While the molecular strain accommo- 
dation mechanism is yet to be established, it manifests as a decrease in the effective 
radius of the molecule and the molecule overwinds initially. Thereafter, the balance 
tilts in favour of the chirality dependent response such that the coupling changes 
its sign and then the molecule begins to unwind. 
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4. Discussion and Conclusions 

The SWCNT and B-DNA response both highlight the link between the microscopic 
strain accommodation behaviour and the macroscopic coupling. In the SWCNT 
case, we have clear evidence that the interplay between chirality and geometry 
is modified by the microscopic energetics, a general design principle for tailoring 
the extent and sign of the coupled response. This could be achieved by chang- 
ing the chemistry of the nanofilamcnt, for example, and has direct relevance for 
tailoring the response of helical nanowires (Kondo & Takayanagi 2000; Tosatti & 
Prestipino 2000), coiled nanobelts (Pan et al. 2001), and nanotube ropes (Liang & 
Upmanyu 2005). In each case, the extent of the response can modify the electrome- 
chanical couplings that can be easily exploited in NEMS and nanoelectronics. Our 
analysis on B-DNA shows that the same design principles are also important in 
scmiflcxiblc, soft chiral filaments. The change in radius of the stretched filament is 
likely to be controlled by the stiffness of the sugar phosphate backbone relative to 
that of the base-pair interactions across the strands. Similar energetic balance may 
also be responsible for the twist-stretch coupling recently extracted in chromatin 
strands (Bancaud et al. 2006), although structural transitions may play a more 
important role there. Identification of such energetics is critical for understanding 
the coupled response in fibrillar assemblies such as coiled-coils and twisted protein 
aggregates. 

An important consequence of the non-linearities is that they also modify the 
effective rigidity and electromagnetic character of these elastic filaments. The ex- 
tensional rigidity of the (8,4) SWCNT develops an additional dependence on the 
induced twist, 

-^ = k e + 2k 12 <p eq . (4.1) 

As a result, the nanotube softens rapidly as it is stretched and is the stiffest when 
compressed to the critical strain. A 2% stretch induces a 7% softening. The effective 
torsional rigidity <9 2 W/<90 2 changes primarily due to the twist induced change in 
chirality. This effect is substantial - a 2% axial strain leads to a 20% change in 
the rigidity (see figure 5A). The reduction the stretch rigidity is also seen in B- 
DNA, where the radius dependence of the rigidities (k e oc R 2 and k<f, cx R 4 ) is 
additionally modified by that of the coupling constant fen (figure 5B), an effect that 
further amplifies the unusually large ratio of the twist to stretching/bending rigidity 
at high strains. Similar softening may also be responsible for the large stretching 
strains that can be absorbed by twisted fibrin polymers (Liu et al. 2006), a crucial 
requirement for the overall viscoclastic properties of blood clots. 

The electromagnetic character of these nanofilaments must depend on the sign 
and extent of these couplings. It is well-known that the SWCNT band gap can be 
tuned by uniaxial and torsional deformation (Minot et al. 2003). Close to the Fermi 
level, the change in band-gap is the highest for achiral nanotubes. The presence of 
this natural coupling in chiral SWCNTs, ignored in these analyses, should increase 
the efficiency of strain engineered band gaps (MU, in preparation). While there is 
considerable debate on the electronic properties of DNA (Dekker & Ratner 2001), 
it is clear that the the inherent coupling must be accounted for in understanding 
the interplay between their flexibility and electronic transport. 
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Figure 5. The influence of twist-stretch coupling on effective rigidities. Strain 
dependence of stretch (left) and torsional (right) rigidities for (A) the chiral SWCNT in 
figure 1 and (B) B-DNA. 

There are several similar examples where Nature has dynamically engineered 
such mechanical couplings. Spiral growth of filamentous algae and fungi (e.g. Phy- 
comyces blakesleeanus) is associated with a sinistral to dextral reversal (Castle 1942; 
Green 1954), a transition in the deformation mode in the fibrillar structure of the 
cell wall in response to turgor pressure (Ortega & Gamow 1974). The sinistral 
growth occurs via reorientation of initially transverse microfibrils along the growth 
direction, while the dextral growth immediately following the sporangium forma- 
tion is dominated by inter-fibril slip (Ortega & Gamow 1974) (figure 6A). The 
microscopic response is essentially what we see in the uniaxially strained SWC- 
NTs. Control studies have shown that the reversal can also be achieved if the cell 
wall is allowed to grow laterally (Yoshida et al. 1980), a combination of the design 
principles underlying the SWCNT and B-DNA response. 

A similar mechanism may also be involved in the propulsion of Spiroplasma, a 
helical wall-less bacteria that infects several plants and insects (Whitcomb 1980). 
The undulations that enable it to swim are due to the formation and propagation of 
kinks, or local reversals in its chirality, along its length (Shaevitz et al. 2005). The 
tapered or bulbous ends of the bacteria can facilitate the nucleation the kinks due to 
motors that generate and then transmit localised axial strains. Another example is 
the cell wall mediated grain spiral decrease and reversal in drying trunks of conifers 
(figure 6B), initiating at the trunk base (Preston 1950; Kubler 1991). Since the 
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Figure 6. Natural examples of engineered reversals in twist-stretch coupling: 

(left) Schematic of spiral growth reversal in late stage (IVa-b) growth of Phycomyces 
blakesleeanus. (inset) Illustration of the rotation reversal due to the transition from fibril 
slip to fibril re-orientation (Ortega & Gamow 1974). (right) Spiral grain variations along 
the stem in a lodgetree pine (Pinus contortd). The white markers are visual guides for 
identifying the change in grain spiral. 

thicker base supports the maximum stress, it likely serves as the origin for changes 
in cell wall microstructure which eventually express as grain spiral reduction and 
reversal. In each case, geometry and the imposed strain effect transitions in the 
interplay between structural chirality and the dominant deformation mode. From 
a teleological perspective, strain engineering such microstructural transitions is a 
powerful instrument allowing dynamic control of the coupling. 
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